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ABSTRACT: In this work we show that the general singlet extension of the MSSM can 
naturally provide a self-interacting singlino dark matter to solve the small cosmological 
scale anomalies (a large Sommerfeld enhancement factor can also be obtained). However, 
we find that the NMSSM (the singlet extension of the MSSM with Z3 symmetry) cannot 
achieve this due to the restricted parameter space. In our analysis we introduce the concept 
of symmetric and antisymmetric viscosity cross sections to deal with the non-relativistic 
Majorana-fermion dark matter scattering. 
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1 Introduction 


As the standard model of the Big Bang cosmology, the ACDM model can account for most 
observations of the Universe. A crucial ingredient of this model is the existence of cold 
dark matter (CDM), which, with a proper cosmological constant, can successfully predict 
the large scale structure of the Universe. However, the predictions on small scale structures 
seem not so successful and some anomalies exist: [1] 1) missing satellites — There should 
be many more dwarf-sized subhalos (satellites) than observed in the DM halo of the Milky 
Way (MW). And the observed galaxy luminosity and Hi-mass functions beyond the MW 
show shallower faint-end slopes than predicted. [2] 2) cusp vs core — It seems to have 
cored inner density profiles in the low surface brightness and dwarf galaxies, this is at odds 
with CDM cusps predicted by simulations [3]. 3) too big to fail — In comparison with the 
densest and most massive satellites found in simulations, the observed brightest satellites 
of the MW attain their maximum circular velocity at a too large radii. [4]. 

There are various ways to solve these small scale problems, such as the nonthermal pro- 
duction of warm dark matter [5] or the baryon feedback in the galaxies to make small halos 
dark [6]. Also, recently the author of [7, 8] proposed another self-interacting Dirac-fermion 
DM scenario with a light mediator (< 100 MeV) to solve these small scale anomalies. 
With a light force carrier, the dark matter scattering cross section could have a non-trivial 
velocity dependence. All of the small scales (the dwarf size, the Milky Way size as well 
as the galaxy cluster size) can have appropriate cross sections, thus leaving enough pa- 
rameter space for the mass of DM, the force carrier and the coupling strength. Besides, 
the authors also showed that the DM self-interactions can be correlated with the effect 
of Sommerfeld enhancement in DM annihilation which is being probed through indirect 
detection experiments. 

This self-interacting DM scenario perfectly explain the anomalies in the simulations 
of small scale structures. Therefore, it is necessary to check if such a scenario can be 
realized in popular new physics theories like low energy supersymmetry (SUSY). In SUSY 
the better known DM candidate is the Majorana-type neutralino, which is composed of 
bino, wino and higgsinos. Apparently, if the neutralino can have self-interactions through a 
light force carrier, such a carrier can not have sizable standard model (SM) interaction due 
to the stringent constraints from both collider and DM detection experiments. Thus, this 
light force carrier should be composed mainly of a singlet with respect to the SM gauge 
groups, which cannot be found in the minimal supersymmetric standard model (MSSM). 

Fortunately, there are various singlet extensions of the MSSM, among which the next- 
to-minimal supersymmetric standard model (NMSSM) seems most attractive [9, 10]. In 
the NMSSM, all the parameters in the superpotential are dimensionless and electroweak 
symmetry breaking is triggered by the TeV-scale soft SUSY breaking terms. The SUSY 
preserving u term in the superpotential of the MSSM is generated by the vacuum expecta- 
tion values (VEV) of a singlet superfield S. It is shown that in the NMSSM a light singlet 
scalar at several GeV can survive the DM detection limits and the collider constraints [11]. 
On the other hand, if we do not impose any discrete symmetry (in the NMSSM it is Z3) 
and allow for all possible interactions of the singlet field, then we have the general singlet 


extension of the MSSM (GMSSM), (more detail can be seen in [12]) which was used to 
explain the PAMELA anomaly [13, 14]. Compared with the NMSSM, the GMSSM has a 
larger parameter space. In the GMSSM, the singlet can form a dark sector in case of a 
very small À. The singlino-like dark matter can annihilate into the light singlet-like scalar, 
which can give the correct DM relic density and a proper Sommerfeld enhancement factor. 
In this model, the singlet scalar can be even lighter than in the NMSSM due to a larger 
parameter space. So it is intriguing to check if such a singlet scalar in the NMSSM or 
GMSSM can serve as the light force carrier mediated in the DM self-interactions, which is 
the aim of this work. 

In this work we focus on the NMSSM and GMSSM to check if the self-interacting DM 
scenario can be realized. In our study we will take into account the constraints from DM 
relic density, the DM direct detection experiments as well as the proper non-relativistic 
scattering cross sections between DM. We organize the content as follows. In Sec. 2, we 
will discuss the general DM interactions. In Sec. 3 and Sec. 4, we will respectively check 
the NMSSM and GMSSM to figure out the possibility of realizing the self-interacting DM 
scenario to solve the small cosmological scale anomalies. Sec. 5 contains our conclusions. 


2 Dark matter interactions 


As mentioned in the introduction, in order to explain both the large scale and small scale 
structures of the Universe, we can introduce the self-interacting DM scenario. The inter- 
actions between DM and SM particles can be summarized as (shown in Fig. 1): 


1. The annihilation to the SM particles (the left diagram of Fig. 1), whose cross section 
at high energy determines the relic density of dark matter and whose cross section at 
low energy is being probed by the indirect detection experiments like PAMELA [19] 
and AMS02 [20]. 


2. The elastic scattering off the SM particles (the middle diagram of Fig. 1), which is 
being probed by various direct detection experiments like CDMS,[15] XENON [16] 
and LUX [17]. 


3. The non-relativistic self-scattering (the right diagram of Fig. 1), where / = 0 in the 
partial wave expansion gives the Sommerfeld enhancement relative to the relativistic 
annihilation while ($25 can account for the anomalies in the small cosmological scales. 


These interactions have subtle correlations among each other. A complete study on DM 
properties needs to combine all these interactions from various experiments. 

We now briefly review the calculation of DM relic density and the DM-nucleon cross 
sections. When the early universe was cooling down, the equilibrium between DM and 
SM particles in the thermal bath can no longer be maintained. The DM will annihilate 
to SM particles until the annihilation rate falls below the expansion rate of the universe. 
Thus the key point in the DM relic density calculation is the annihilation rate of DM. 
In our following calculation, we use the standard method [18] to calculate the relativistic 
annihilation cross section and the degrees of freedom at the freezing out temperature. 
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Figure 1. Dark matter interactions: annihilation to SM particles in the left diagram, scattering 
off quarks in the middle diagram, and self-scattering in the right diagram. 


It is sufficient to consider only the spin independent (SI) elastic cross sections between 
DM (denoted by x) and nucleon (denoted by fp for proton and fn for neutron [18] ) because 
of its high sensitivity in current DM direct detection experiments. These interactions are 
dominantly induced by scalar exchange processes at tree level, as shown in the middle 
diagram of Fig.(1). Note that the vector boson exchange interactions are also possible and 
we concentrate on the scalar interactions only. For moderately light scalar bosons, fp is 
approximated by (similarly for fn) 


Ko w2 i 

f= D, mii tofe 2, mp, (2.1) 
m q 27 m 
q=u,d,s 1 q=c,b,t 4 


where fee denotes the fraction of mp (proton mass) from the light quark q, fr, = 1 — 
i de fe ) is the heavy quark contribution through gluon exchange, and fe is the coefficient 
of the effective scalar operator given by 


fe 2 Coxa, (2.2) 


mg 
with C being the corresponding interaction vertex and mg the mass of the exchanging 
particle. The DM-nucleus scattering rate is then given by 


2 
where my is the DM mass, mr is the mass of target nucleus, and np(nn) is the number of 
protons (neutrons) in the target nucleus. We can see the dependence of SI cross section 
on the mass of exchanging particle o% œ 1/ ms which is very important in our following 
discussions. 

Note that the annihilation cross section in relic density is calculated at high energy. 
To explain the DM indirect detection results (such as the PAMELA result) and the small 
cosmological scale anomalies, the DM scattering at low energies are needed. We will discuss 
their relations in the following. 


2.1 The Sommerfeld enhancement effect 


The Sommerfeld enhancement effect in dark matter annihilation is proposed to explain 
some DM indirect detection results, such as the positron excess observed by PAMELA [19] 


or AMS [20]. The explanation of positron excess requires a very large DM annihilation rate 
which on the other hand can not explain the DM relic abundance if the DM is produced 
thermally in the early universe. The Sommerfeld effect can greatly enhance the annihilation 
rate when the velocity of DM is much smaller than the velocity at freeze-out temperature. 
Note that it is shown that the positron excess can be more naturally explained by the pulsar 
wind and, further, the gamma-rays produced by inverse Compton scattering on interstellar 
radiation field of electrons and positrons produced by dark matter basically exclude the 
dark matter interpretation [21]. 

The Sommerfeld enhancement is a common effect in non-relativistic quantum mechan- 
ics (QM). When DM has self-interactions, it will lead to an effective potential V(r) in the 
non-relativistic limit. Thus the Schrödinger equation of DM particle can be written as 


2m 


Te. k? 
— — Vy, + Vind, = gm de (2.4) 
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where k is the relative momentum of the DM particle. A non-relativistic DM moves and 
annihilates around the origin, namely near r = 0. Therefore, the only effect of the potential 
V(r) is the change of the modulus for the wave-function at the origin comparing to that 
without V(r). Then, the annihilation cross section is enhanced to 


0 = Tosk, (2.5) 


where the Sommerfeld enhancement factor S' is given by 


_ WO? _ 1, yp 
s= aye O (2.6) 


with ~,(0) and po (0) are respectively the perturbed (unperturbed) wave function of DM. 

We can find a solution for the Schrödinger equation in non-relativistic QM. As Ry(r) ~ 
rl as r > 0 and the l > 0 part will not change the annihilating wave function at the origin, 
we need only to calculate the l = 0 part of the radial wave Ryi(r). We can numerically 


solve the Schrödinger equation with the boundary condition 
x(r) = rRe(r) > sin(kr +ô) as r— oo. (2.7) 


In case of a Yukawa-type effective potential induced by scalar exchange 


V(r) = -ay Z= ; (2.8) 


r 


with ay = |Cgy,|?/4m, the Sommerfeld enhancement factor is given by 


(2.9) 


Though there is no analytical solution for the Schrödinger equation with the Yukawa- 
type potential, there are three distinguishable regions for the Sommerfeld enhancement, 
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Figure 2. The Sommerfeld enhancement factor at different velocities with the coupling strength 
ay = 0.01. Here v = 10km/s, 200km/s, 1000km/s correspond to the characteristic speed of dwarf 
halos, the Milky Way and clusters, respectively. 


depending on the value of m,/mg. If the mediator scalar mass is comparable to the 
DM mass, Sz is negligible at all scales. However, if the mediator scalar mass is much 
smaller than the DM mass, the enhancement factor becomes independent of m,/mg and 
corresponds to the Coulomb limit. In this limit, the Sommerfeld enhancement factor is 
essentially given by Sy ~ ma,/v. In the resonance regions where my/mg ~ 7?n?/6ay 
with n = 1,2,3,---, the DM annihilation cross section can be enhanced. The enhancement 
factor is approximately given by Sk ~ Tay Mg / (6m, v7), which is very sensitive to the DM 
velocity. In our numerical calculation we reproduced the results of [23], as shown in Fig. 
2. This figure shows the Sommerfeld enhancement factor at velocity 10, 200, 1000 km/s, 
corresponding to the characteristic velocities of the halo of the dwarf, the Milky Way and 
the cluster (these are the three small cosmological scales at which the ACDM model seems 
not work well). 


2.2 The transfer and viscosity cross sections of self-interacting dark matter 


As pointed in [23, 24], in order to solve the small scale simulation anomalies, self-interaction 
between the DM is necessary. In the non-relativistic limit, the scattering between DM can 
be described by QM. The most recent simulations have shown that ¢/m, ~ 0.1—10 cm?/g 
on dwarf scales (the characteristic velocity is 10 km/s) is sufficient to solve the core-vs-cusp 
and too-big-to-fail problems, while the Milky Way (the characteristic velocity is 200 km/s) 
and cluster scales (the characteristic velocity is 1000 km/s) require g/m, ~ 0.1—1 cm?/g. 
It appears that all the data may be accounted for with a constant scattering cross section 
around o/m, ~ 0.5 cm?/g. On the other hand, the self-interacting DM models generically 
predict a velocity-dependent scattering cross section over a wide range of parameter space. 

The numerical input for the simulation of small scales is the differential cross sec- 
tion do/dQ as a function of the DM relative velocity v. The standard cross section 
o = f dQ(do/dQ) receives a strong enhancement in the forward-scattering limit (cos 0 — 1), 


which does not change the DM particle trajectories. Thus two additional cross sections are 
defined to parameterize transport [25], the transfer cross section øp and the viscosity (or 
conductivity) cross section ov: 


or = | a@ (1 — cose) F ov = [ao sin? Z (2.10) 
The transfer cross section is weighted by (1 — cos@) which is the fractional longitudinal 
momentum transfer while the viscosity cross section is weighted by sin? which is the 
energy transfer in the transverse direction. The transfer cross section was used in the 
DM literature to regulate the forward-scattering divergence in case of Dirac-fermion DM 
candidate, while the viscosity cross section was used in case of Majorana-fermion DM 
candidate. This is because both forward and backward scatterings diverge, corresponding 
to poles in the t- and u-channel diagrams for the identical DM candidate. 

In [23] the authors considered a Dirac-fermion DM and calculated or. Within the 
resonance region, no analytic formula exists for op and it must be computed by solving 
the Schrodinger equation directly with partial wave expansion method. The scattering 
amplitude is given by 


f@)= DDE + 1)et® P;(cos 0) sin dg . (2.11) 


The differential scattering cross section is given by 


do Lj i 2 

qa = a So (26 + 1)et™ P;(cos 0) sin de), (2.12) 
1=0 

where 6, is the phase shift for a partial wave £. In terms of the phase shifts, the transfer 

cross section is given by 


ork? 2 . 

aie NO (2+ 1) sin? (Se41 — 8e) . (2.13) 
é=0 

To obtain 6g, one must solve the Schrödinger equation to calculate the radial wave function 

Re(r) for the reduced two-DM particle system 
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The Schrödinger equation can be recast into the form 


@  , +I l op 
(a+ = ge ae ) xela) = 0 (2.15) 


with new variables 


Verh, rt=aymyr, a=;—, b= (2.16) 
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Figure 3. The birds-eye view of transfer and viscosity cross sections or, oy s, ov 4 in the parameter 
space (a, b). 


When the DM candidate is a Majorana fermion, the amplitude and Schrödinger equa- 
tion are the same as in Eq.(2.11) and Eq.(2.15). However, the total wave function of the 
spin-1/2 fermionic DM must be antisymmetric with respect to the interchange of two iden- 
tical particles. Then the spatial wave function should be symmetric when the total spin is 
0 (singlet) while the spatial wave function should be antisymmetric when the total spin is 
1 (triplet). The viscosity cross section should be defined with two variables: 


2 
d 1 > l 
E SAO- =z] DO (22+ 1)(exp(2idy) — 1) Pe(cos 8)| (2.17) 
¢(EVEN number) 
2 
d 1 > , 
A =|f@-Ff-OP =| JO (22+ (exp(2ið) — 1)P:(cos6)| (2.18) 
£(ODD number) 


Using the orthogonality relation for the Legendre polynomials, we can obtain 


k? = 
avs = XO 4 sin? (ep2 — ôe) (E + 1)(€ + 2)/ (22 +3), (2.19) 
i L(EVEN number) 
k? °° 
wa = So Asin? (Sep — 5¢)(6+1)(¢ + 2)/ (26 + 3). (2.20) 


(ODD number) 


From the expressions of the transfer and viscosity cross sections, we can see that both op 
and oy will converge to a static value as the phase shift 6¢ approaching to a same value 
when the partial wave @ grows up. We adopt the numerical method in [23] to calculate all 
the cross sections. The results are shown in Fig. 3 in which we have the same definitions 
of input variables. The left plot is the reproduced results, same as in [23]. We can see that 
or and oys are almost same while oy 4 shows differences in some region (ab < 0.1). The 
ab > 0.5 and b > 1 is the classical regime in which all the cross sections can be very large. 
There are also the resonance regions (b = 1 and ab < 0.5) in which all the cross sections 


exhibit patterns of resonance, making the cross section much complicated. Self-interacting 
DM scenario used in [23] and in this paper works around the resonance regions. In the 
region ab < 0.1, oy, is much smaller than the other two cross sections. The reason is that 
the summation of the phase shift in oy 4 begins at @ = 1 as shown in Eq. (2.20), and in 
this region the dominant phase shift is 69 and other phase shifts are much smaller. That 
is why we always calculate 69 for partial wave expansion in QM. In our following analysis, 
we assume that the DM scatters with random orientations, thus the triplet is three times 
as likely as the singlet and the average cross section will be 

oy = aes F -a ; (2.21) 

4 4 

This is simple but sufficient to estimate the viscosity cross section for Majorana DM. 


3 Can self-interacting dark matter scenario be realized in the NMSSM 
? 


SUSY can not only give a solution to the hierarchy problem but also provide a good dark 
matter candidate and realize gauge coupling unification. Among the SUSY models, the 
MSSM has been intensively studied. This model, however, has the little hierarchy problem 
since the newly discovered 125 GeV Higgs boson requires a heavy stop or a large trilinear 
coupling A;. Besides, the MSSM also suffers from the -problem [26]. It is remarkable that 
both the little hierarchy problem and the u-problem can be solved in the NMSSM (9, 10], in 
which an additional gauge singlet ŞS is introduced. In this model the p-term is dynamically 
generated through the coupling SH„Ha after S develops an electroweak scale VEV, while 
the little hierarchy problem is solved through an additional tree-level contribution to the 
Higgs mass. With the additional singlet, it might be possible for the NMSSM to give DM 
a proper non-relativistic cross section by tuning the singlet mediator. In the following we 
check this possibility. 


3.1 Dark matter and Higgs bosons in the NMSSM 
In the NMSSM the relevant superpotential containing S is given by 

AGH ys Hy t =o (3.1) 
where H,, and Hy are the Higgs doublet superfields, and A and « are dimensionless pa- 
rameters. Note that there is no explicit w-term and an effective u-parameter is generated 


when the scalar component (S) of S develops a VEV s: pe ff = As. The corresponding soft 
SUSY breaking terms are given by 


A ASH, - Ha + fens? + h.c.. (3.2) 
So the scalar Higgs potential is given by 
Vr = |AHa: Hy — 6S?|? + JAS)? (Hal? + |Hul?) (3.3) 
Vo = È (Ma P1 al? — IHa BaP) +222 Qaa- P), (3.4) 
Veors = mil Hal? + mi] Hl? + m2|SP? — (AASHa- Hu + AnS? + he.) , (3.5) 


where gı and go are the coupling constant of Uy(1) and SUz(2), respectively. Assuming 


H, iH, H H S LS 

uR t? ul dRt+? at S=s4 RTI (3.6) 
v2 v2 v2 
with hu, ha, s being the corresponding VEVs and using the minimization conditions, one 
can obtain a 3 x 3 CP-even Higgs matrix Mp, a 3 x 3 CP-odd Higgs matrix Ma and a 


2 x 2 charged Higgs matrix Me. Note that there are three Goldstone bosons in Ma and 


H? = hy + Hy = ha + 


M, which imply that after diagonalization these two matrices must have 0 eigenvalues. 

From the superpotential in Eq. (3.1), we can see that the interactions between the 
singlet and the SM sector are controlled by the parameter A. If a light singlet Higgs exists, 
the constraints from both collider and DM detections can be satisfied only if A is small 
enough. Then the singlet Higgs will be a dark sector. The spectrum of the NMSSM has 
been widely studied in the literature [10], so we only present our conventions and list two 
necessary terms for our analysis, concentrating on the dark singlet sector: 


e The CP-even Higgs mass matrix Mp in the interaction basis $°¢"* = (Hur, Har, SR) 
can be diagonalized by an orthogonal matrix S;; to obtain 3 CP-even mass eigenstates 
(ordered in mass) h; = Si, Share with the corresponding masses denoted by mp,. The 
elements involve the singlet component are 


uh 
Mnh,33 = TRA can KS(A, + 4K), (3.7) 
Mnh,13 = 2A Hef hu = Ahal Aa + 2ks), 
Mh,23 = 2A Lett ha = Ahul Ay + 2ks). (3.9) 


e Neutralino matrix Mo is composed of the U(1)y gaugino \;, the neutral SU (2) 
gaugino A2, the singlino ~, and the neutral higgsinos Waa: In the basis y? = 
(—iA1, —ià2, YL, Y, Ys) one can rewrite 


1 
l= (Y?) TMo (yY?) + h.c. (3.10) 
where 3 E 
gıhu _ Jiha 
M, 0 v2 5 0 
g2hu g2hd 
o Me Ee o 
Mo =| 27 ec 0 per —Ahg |- (3.11) 
_gha gzħa _ 0 a 
V2 J2 Leff u 
0 0 —Ahg —Ahy 2ks 
After diagolization, one obtains 5 eigenstates (ordered in mass) x? = ij 9, among 


which the lightest one x? is usually assumed to be the lightest SUSY particle and 
serve as a good DM candidate. 


From the spectrum we can see that if A approaches zero, one can tune the masses of the 
singlet-dominant scalars and neutralinos to arbitrary values by varying the parameters 
k, Ag. 


—10— 


DM (the lightest neutralino) in the NMSSM has three-type components: gaugino, 
higgsino and singlino. Assuming the gaugino unification relation Mj/M, ~ 2, we have 
three possibilities for the DM: 


e Bino-dominant DM. As shown in [27], under current collider and DM relic density 
constraints, the SI cross section can exclude a large part of parameter space, leaving 
only a bino-dominant DM candidate below TeV. 


e Higgsino-dominant DM. As pointed in [28], the higgsino-dominant DM candidate 
around 1.1 TeV can satisfy all the constraints, including the relic density and current 
DM direct detections. 


e Singlet-dominant DM. In order to explain the observation of CoGeNT [29], the analy- 
sis in [30] showed that in the Peccei-Quinn limit there can exist three light singlet-like 
particles (0.1-10 GeV): a scalar, a pseudoscalar and a singlino-like DM candidate. 
For a certain parameter window, through annihilation into the light pseudoscalar 
the singlino DM can give the correct relic density, and through exchanging the light 
scalar in scattering off the nucleon a large cross section suggested by CoGeNT and 
DAMA/LIBRA [31] can be attained. 


3.2 Allowed parameter space and dark matter self scattering in the NMSSM 


We use the package NMSSM Tools [32] to study all the scenarios numerically. As the package 
requires the mass of light Higgs to be heavier than 1 GeV, we modify the package so that the 
light Higgs mass can be arbitrarily low and also we include our own codes in the calculation 
of SI cross section. Since it is nontrivial to obtain a light singlet-dominant Higgs, we 
randomly scan the parameter space under the condition mp, < 5GeV and mp, > 120GeV. 
In our scan we consider constraints from DM relic density and DM-nucleon scattering cross 
section. Since the SM Higgs (h2) is much heavier than the singlet-like Higgs (h1), we only 
need to slightly tune the parameter A, to get a light Higgs with mass 100 MeV. Finally, we 
use a light scalar with mp, = 100MeV for the calculation of DM-nucleon SI cross section 


SI. Note that the singlet sector will be dark with respect to the SM sector when A is 


o 
sufficiently small, and then the SM sector is very insensitive to the tuning in the dark 
sector. 

The numerical results of the three possibilities are shown in Fig. 4, in which the SI 
cross section of DM-nucleon scattering is shown and current constraints of XENON and 
LUX are displayed. We can see that the DM direct detection constraints exclude most 
of the samples. However, there are still some surviving parameter space, especially for 
bino-dominant and higgsino-dominant scenarios. Future XENONIT will further cover the 
surviving parameter space. 

Now we present some details related to our scan for different scenarios: 


e For the singlino-dominant DM scenario, we define the parameter £ = Au/mz, € = 
A)/ptan 8 — 1 and scan the parameter space as in [30]: 


2 < tan 8 < 50, 0.05 < A < 0.2, 0.0005 < « < 0.05,-0.1 < £ < 0.1, 
ewe',|A,x| < 500GeV, |u| < 1TeV. (3.12) 
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Figure 4. The spin-independent cross section of DM scattering off the proton versus the DM mass 
in different scenarios in the NMSSM. [O denotes singlino-dominant DM scenario, most of which are 
excluded by the DM direct detection limits; A denotes bino-dominant DM scenario, among which 
the red ones are excluded while the green ones survive; o (around 1 TeV) denotes the higgsino- 
dominant DM scenario, among which the pink ones are excluded while the green ones survive. All 
the samples satisfy the DM relic density constraints. 


The sfermion mass is set to be 6 TeV so that we can easily get a 125 GeV SM Higgs. 
As noted in [30], the DM annihilates to SM particles mainly through the resonance of 
the singlet pseudo-scalar a; (the Feynman diagram is shown in Fig. 5). The result of 
the scan is shown in the left panel of Fig. 6. As a consequence of resonance, the singlet 
pseudo-scalar mass is almost twice of the DM mass. This result is consistent with 
[30] and also implies that the coupling of a, f f should not be completely neglected 
and thus the parameter A can not be too small. On the other hand, as shown in sec.2, 
the SI cross section has the inverse quartic power dependence on the mediated scalar 
mass (here it is the singlet Higgs h1). Though the coupling of hı f f is suppressed by 
a small A, the cross section will be enhanced greatly when the mass of hı is below 
GeV. Therefore, this scenario can hardly satisfy the direct detection limits on the SI 


cross section. 
e For the bino-dominant DM scenario, we scan the parameter space 
1075 < JA, |x| <1, 2<tanB<50, |py|,|Ay|<1TeV, |M2| < 500GeV , (3.13) 


and set A, ~ —4K/X to get a light singlet Higgs (see Eq. (3.8)). We also keep the 
gaugino unification assumption Mı = M2/2, M3 = 3Mp2 and set the soft parameters 
to 6 TeV. 
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Figure 5. Feynman Diagram of DM annihilation to SM fermions through a pseudo-scalar a. 
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Figure 6. The scatter plots of parameter space allowed by DM relic density shown in the plane 
of Ma, versus the DM mass. The left panel is the singlino-dominant DM scenario, where most of 
the samples are excluded points by direct detection limits of LUX-2013. The right panel is the 
bino-dominant DM scenario, where the red ones are excluded points by direct detection limits of 
LUX-2013 while the green ones survived. 


Since bino does not couple to singlet Higgs directly, the coupling between bino- 
dominant DM and singlet Higgs arises from its higgsino and singlino components. 
Therefore, the SI cross section of bino-dominant DM can be much less than the LUX 


limits, as shown in Fig. 4. 


The right panel of Fig. 6 shows the samples that satisfy the DM relic density in the 
bino-dominant DM scenario. The samples can be divided into two regions: one is in 
the diagonal region(m, ~ 2m,) and the other is a vertical belt around my ~ 50GeV. 
The first region is excluded by the direct detection limits of LUX (the reason is similar 
to the singlino-dominant DM scenario). In the second region, the DM annihilates to 
SM fermions through the Z-boson resonance, independent of singlet coupling A. The 
coupling A is constrained only by the SI cross section limits and can only take a very 
small value, namely \<0.005. This can be seen from the left panel of Fig. 7. 
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Figure 7. Same as Fig.4, but showing the spin-independent cross section versus À in the NMSSM. 


e For the higgsino-dominant DM scenario, we scan the parameter space 


107° < |Al, |x| <1, 2< tang <50, |A| <1Tev, 
Jul <1.5TeV, Ay ~ —4ku/à, |Mo| = 10TeV . (3.14) 


As above, the gaugino unification assumption is used and the soft parameters are set 
to be 6 TeV. In this scenario, the DM mass need to be around 1100 GeV to satisfy 
the relic density, and thus no resonance through Z-boson or pseudo-scalar happens 
in the DM annihilation. 


The singlet Higgs can couple directly to the higgsino through the AS Yupa term in 
the Lagrangian, so the SI cross section limits only constrain A. From the right panel 
of Fig. 7 we can see that A<0.001 for the survived samples. 


As a brief summary, the singlino-dominant DM scenario is almost excluded because 
of the A-value correlation between relic density and DM-nucleon SI cross section (the relic 
density needs a not-too-small which is too large for the DM-nucleon SI cross section). 
On the other hand, some parameter space still survived the SI cross section limits in 
the bino-dominant and higgsino-dominant DM scenarios, in which the \-value correlation 
between relic density and DM-nucleon SI cross section is relaxed and A is constrained by 
the DM-nucleon SI cross section limits to take a very small value. 

In order to figure out if the NMSSM can realize the DM self-interaction scenario to 
explain the small cosmological scale anomalies, we now calculate the DM self scattering 
cross section. For each allowed sample in Fig. 7, we use three velocities corresponding to 
the three small cosmological scales (dwarf, Milky Way and cluster) to calculate the DM self 
scattering cross section. Then for each sample we have three values of DM self scattering 
cross section and we keep the maximal one to display it in Fig. 8. We see that due to a 
very small À value, oy /m, is too small to explain the small cosmological scale anomalies 
(its value must be in 0.1 — 10 cm?/g to explain dwarf scale anomaly and in 0.1 — 1 cm?/g 
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to explain MW and cluster anomalies). For such a weak scattering, we can regard the DM 
as a collisionless particle. So we conclude that NMSSM can not realize the self-interacting 
DM scenario to solve the small cosmological scale anomalies. 
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Figure 8. Same as Fig.7, but showing the non-relativistic viscosity cross section between DM in 
the NMSSM. For each allowed sample in Fig. 7, we use three velocities corresponding to the three 
small cosmological scales (dwarf, Milky Way and cluster) to calculate the cross section and here we 
plot the maximal one. 


4 Can self-interacting dark matter scenario be realized in the GMSSM 
? 


The main reason for the failure of the NMSSM in solving the small cosmological scale 
anomalies is a too small A (constrained by DM-nucleon SI cross section limits), which 
determines the coupling strength of hıxx in DM self-interaction (Note that when the DM 
is singlino-dominant, the coupling hıxx can be proportional to «x. However, a singlino- 
dominant DM is obtained under the condition k < A [30], which means a very small x). 
In the NMSSM the ASH, - Hq term in the superpotential is the origin of the p-term. 
In the GMSSM, however, the Z3 discrete symmetry is not imposed and the u term can 
exist in the superpotential, together with the ASH, - Hg term (several other terms of the 
singlet superfield can also exist in the superpotential). Then the singlet sector can be a 
completely dark sector in case of a nearly vanishing A. Unlike the NMSSM, the dark Higgs 
sector (including a singlino-dominant DM) can be easily realized in the GMSSM, which 
does not need the condition x < A [30]. This means that a singlino-dominant DM can be 
obtained with a sizable « and in this case the coupling hıxx in DM self-interaction, which 
is proportional to «, can be large. So it should be easy to realize the self-interacting DM 
scenario in GMSSM. 
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As shown in [13], in the GMSSM a singlino-dominant dark matter annihilates to a 
singlet-dominant Higgs, which can give the correct DM relic density, and also the DM- 
nucleon SI cross section limits from the direct detections can be easily satisfied. In addi- 
tion, an appropriate Sommerfeld enhancement can explain the positron excess observed by 
PAMELA. The non-relativistic DM scattering and Sommerfeld enhancement are similar 
to what discussed in [7] except that there can mediated an additional pseudo-scalar in the 
GMSSM. 


4.1 Dark matter and Higgs bosons in the GMSSM 
The superpotential of the GMSSM is 


we A ee x l. å 1 A 
W = uH,- Ha + ASH,- Hat+nS + a TAn m ; (4.1) 


which involves the parameters u, A, N, Hs, kK. We set A ~ 0 so that the singlet sector 
almost decouple from the SM sector. In the following discussion, we will concentrate on 
the singlet sector. 

The soft SUSY breaking terms take the form 


1 1 
= Leo, = m2|s|? + (cans + 5 Bigs” + 3 KA, So + he.) . (4.2) 


The main difference between the NMSSM and GMSSM is reflected in their Higgs sectors 
which contain different singlet Higgs mass matrices and self-interactions. The difference 
mainly comes from the F-term F;: 


Vr, = |Fsl? = |n + s8 + «S|? 
= | KS? +n? + p2|S|? + (nusS + KNS? + Kus S? S* + h.c.). (4.3) 


Since 7? is a constant, the u2|S|?, nus 9, «nS? terms can be absorbed by the redefinition of 
the soft SUSY breaking parameters m2|S|?, C,nS, 5 BsusS?. 
The singlet Higgs potential is 


V = Vr + Voott 
1 1 
m2|$|? + (cons 5 B,ji,57 + 3 An SÈ + Kus S? S* + he.) . (4.4) 


The singlet chiral supermultiplet contains a complex scalar and a Majorana fermion y. 
After the scalar component getting a VEV, we can get one CP-even Higgs h and one 
CP-odd Higgs a. The mass spectrum and the relevant Feynman rules are 


My = 2ks + pus, (4.5) 
m? = ks(4ks + A, + 3u5) — Cyn/v, (4.6) 
m2 = —2BsHs — 68(3A, + us) — Cyn/s, (4.7) 
Vann = —4ik(6Ks + Ay + 35) = —4in(3m, + Ax), (4.8) 
Vhaa = —4ik(2Ks — A, + us) = —4in(m, — Ax), (4.9) 
Vhyx = —4tk, (4.10) 
Vaxy = —4677°. (4.11) 
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We do not give the Feynman rules for the four-scalar vertex for they are irrelevant to 
our analysis. From the spectrum and Feynman rules above, we can see that the mass 
parameters of the three particles (h,a, x) can be set freely because the relevant input soft 
parameters can take arbitrary values. This makes the following calculation much easier 
and we can choose the input parameters easily without fine-tuning. 

Then we have only five input parameters in our calculation, namely my, Mp, Ma, K, Ap- 
It is useful to note that the three CP-even Higgs vertex Vpnn can vanish if we set A, = 
—3my,. 


Figure 9. Feynman diagram of singlino DM annihilation to scalars in the GMSSM. The final state 
of ¢ ¢ can be Ah, ha and aa. 


The calculation of non-relativistic DM scattering cross section is already given in Sec. 
2 with the coupling strength of DM to mediator (a, = 4K*/m). Here we just give some 
details relevant to the calculation of DM relic density. The Feynman diagrams of the DM 
annihilation are shown in Fig. 9. The general form of the annihilation cross section is given 
by [33] 7 
ei 
48rsS 
where S$ is the symmetry factor, A('S9), A( Po), A Pı), A(? P2) are the contributions 
from different spin states of DM and (y is given by 


By = y/1 = 2(m3 + mp) + (m3 + 3/2/92 , (4.13) 


with X,Y being the final states. The amplitudes from different final states are given by 


Ov 


ACSO + 5 (ACRE + ACP IP) +ACPAP|, (412) 


1. xx> hh: 
R(3my + Ax) 1+R(m) 457 
A(?P9) = 1 a EE 
(Py) = 16v Gun" | R(mp)2 + iGh PB 3P2|’ ae 
A(??P2) = —(128/V3)un?83/P2 . (4.15) 
2. xx 7 aa 
R(m, — Ax) 1—R(m,) 46 
ACH = 1 2 | GE ee 
(Po) = 16v6un? | R(mp)2 + iGh P, 3 P2| ee) 
A(? P2) = —(128/V3)un? 8} / P? . (4.17) 


= 7 = 


3. xX > ha: 


R(my) 1 1 Bs 
AV 2K? X |14 | =- — H Ml 
ave B |; 2P, 3P? 
2 2 2 2 1 1 BF 
um Ao 14 See a l A 
16V 2k (R(ma)? — R(mn)*) | v 8 2P, * BP (4.18) 
AČ Pi) = 64uK?B3/P2 . (4.19) 
In the above formulas 
m 1 Irem 
R(mx) = —, Py =1+ R(m,)? — =(R(mx)? + R(my)}?), Gn = 5%. (4.20) 
My 2 sie 


Given the annihilation amplitudes, the relic density Qh? can be calculated in the standard 
way. From the expression of the DM relic density [23] 


3 Ta? A 
ov = 2 Xy, |1 — h ' (4.21) 
8 m? m? 


we can see that the ratio a,/m, is almost a constant when mp < my in order to give the 
correct DM relic density. 


4.2 Allowed parameter space and dark matter self scattering in the GMSSM 


We first change our Feynman rules to those in [23] and check if we can reproduce the 
results in [23]. Our results are shown in Fig. 10, which agree with [23]. Fig. 10 shows the 
transfer cross section or allowed by solving the small cosmological scale anomalies (dwarf, 
Milky Way and cluster) with their corresponding characteristic velocities. Here all the 
samples satisfy the DM relic density given by PLANCK [34]. We can see that o7/my can 
satisfy simultaneously the requirements of dwarf, Milky Way and galaxy cluster scales to 
solve the small scale anomalies. On the other hand, as can be seen from the figure, the 
requirement to solve such small scale structure anomalies can also give strong constraints 
on the parameter space of self-interacting DM scenario. This means that subtle relations 
should be satisfied between mp, my and ay. 

Then we perform a numerical calculation in the GMSSM. Since in the GMSSM the 
singlino DM is of Majorana type, we need to use viscosity cross section oy. In Fig. 11 we 
show the GMSSM parameter space allowed by DM relic density plus all three small scale 
structures. The corresponding parameter space for the DM self-interaction model [23] is 
also displayed for comparison. We can see that compared with the DM self-interaction 
model [23], the GMSSM has a larger parameter space to solve the anomalies of all three 
small scales (the mass of DM can be heavier than 10 GeV and the coupling strength a, can 
be at order 107°). The reason is that in the DM self-interaction model [23] DM can only 
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Figure 10. Scatter plots of the parameter space in the DM self-interaction model [23] allowed by 
DM relic density plus dwarf scale structures (blue region), Milky Way size structures (red region) 
or cluster size structures (green region). We obtained these results by changing our Feynman rules 
to those in [23]. 
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Figure 11. The scatter plots of the parameter space allowed by DM relic density and all three 
small scale structures (dwarf, Milky Way, cluster). The red ones are for the DM self-interaction 
model [23] while the green ones are for the GMSSM. 


annihilate into hh via t-channel and u-channel while in the GMSSM DM can annihilate 


= {9 = 


10° l 
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10 10° 


Figure 12. Same as Fig.10, but showing the Sommerfeld enhancement factor. The left frame is 
for the DM self-interaction model [23] while the right frame is for the GMSSM. 


into hh, ha and aa via t-channel, u-channel and s-channel, as shown in Fig.9. 

Finally, in Fig. 12 we show the Sommerfeld enhancement in the allowed parameter 
space. We see that compared with the DM self-interaction model [23], the GMSSM allows 
for a larger m,/mp, and thus can give a larger Sommerfeld enhancement factor. The 
ongoing indirect detections of DM can probe the DM annihilation rate (shed light on the 
Sommerfeld enhancement factor) and thus help to distinguish different DM self-interaction 
models. 


5 Conclusions 


In this paper we studed the possibility of the DM self-interaction to solve the small cosmo- 
logical scale anomalies in the singlet extensions of the MSSM. We first checked the NMSSM 
and found that the correlation between the DM annihilation rate and DM-nucleon SI cross 
section strongly constrains this model so that it cannot realize the DM self-interacting sce- 
nario. For the GMSSM, the parameter space was found to be large enough to realize the 
DM self-interacting scenario and at the same time can give a large Sommerfeld enhance- 
ment factor. Also, we found that for Majorana-fermion DM, we must use viscosity cross 


sections (oys and oya) in DM simulations. 
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